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Abstract. Explicit solutions to the Riemann-Hilbert problem will be found realising 
some irreducible non-rigid local systems. The relation to isomonodromy and the sixth 
Painleve equation will be described. 



1. Introduction 

Unfortunately, to say that a particular Riemann-Hilbert problem is "solvable", one 
usually means that there exists a solution rather than that one is actually able to solve 
the problem explicitly. 

In this article we will confront the problem of explicitly solving the Riemann-Hilbert 
problem directly, for irreducible representations (so we already know the problem is "solv- 
able"). We will describe how one soon becomes embroiled in isomonodromic deformation 
equations, from which it is easy to see the difficulty: in the simplest non-trivial case the 
isomonodromy equations reduce to the sixth Painleve equation Pyi and one knows that 
generic solutions of Pyi cannot be written explicitly in terms of classical special functions. 

However there are some explicit solutions to Pyi and our aim will be to write down some 
new solutions controlling isomonodromic deformations of non-rigid rank two Fuchsian 
systems on the four-punctured sphere. The cases we will study here will have monodromy 
group equal to either the binary tetrahedral or octahedral group (the icosahedral case 
having been studied in 0), or to one of the triangle groups A237 or A238. 

Previously a tetrahedral and an octahedral solution of Pyi have been constructed by 
Hitchin ^31 and (up to equivalence) independently by Dubrovin jH]. Moreover with hind- 
sight we see there are three other such solutions in the work of Andreev and Kitaev PJ QBl • 
Here we will classify all such solutions and find an explicit solution in each of the new 
cases that appear. 

Amongst the solutions which look to be new (i.e. to the best of the author's knowledge 
have not previously appeared) there are five octahedral solutions including one of genus 
one, and two 18 branch genus one solutions with monodromy group A237. The largest 
octahedral solution has sixteen branches which is (currently) the largest known genus 
zero solution (those with more branches in j3] having higher genus) and we will show it 
is equivalent to a solution with monodromy group A238- 

The results of sections 0] and IH will be of particular interest to people interested in 
constructing linear differential equations with algebraic solutions (cf. e.g. jTH El EZl E! ) • 
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Indeed tables 1 and 3 may be interpreted as the analogue for rank two Fuchsian systems 
with four poles on P^, of the tetrahedral and octahedral parts of Schwarz's famous list 
P3] of hypergeometric equations with algebraic bases of solutions. 

2. From Riemann-Hilbert to Painleve 

Consider a logarithmic connection V on the trivial rank n complex vector bundle over 
the Riemann sphere with singularities at points Choosing a coordinate z 

on the sphere (in which a^n = oo say), this amounts to giving the Fuchsian system of 
differential equations V d/dz which will have the form: 

, m—l . 

(1) ;r-^^^)' ^W = E^ 

for complex nxn matrices Ai. The original Riemann-Hilbert map is the map which takes 
such a Fuchsian system to its monodromy data: restricting V to the punctured sphere 

P* :=pi\{ai,...,a„} 

yields a nonsingular holomorphic connection and taking its monodromy yields a repre- 
sentation 

p G Hom(7ri(P*),G) 

where G = GL„(C). The Riemann-Hilbert problem is the following: given ai, . . . , and 
p can we find such a connection V with monodromy equal to p? 

Upon choosing simple loops 7^ in P* around generating 7ri(P*) and such that 7mO- ■ -071 
is contractible one sees that for each m-tuple of points a = (ai, . . . , a^) the Riemann- 
Hilbert map amounts to a map between the following spaces: 

(2) {{A,,...,AJ I 5^A, = 0} ^ {(Mi,...,Mj I M^---Mi = l} 

where Mi = p{'~fi) G G. The Riemann-Hilbert problem then becomes: given a point 
M = (Ml, . . . , Mm) on the RHS of are there matrices A = {Ai, . . . A^) with ^ A^, = 
on the LHS such that RHa(A) = M? 

Remark 1. So far we have ignored the questions of choosing a basepoint for 7ri(P*) and 
the choice of basis of the fibre at the basepoint. However it is immediate that if we have 
a solution RHa(A) = M (defined with respect to some choice of basepoint /basis) then 
conjugating the matrices Ai by some constant matrix g E G corresponds to conjugating 
the monodromy matrices Mj as well. Thus the Riemann-Hilbert problem is independent 
of the choice of basepoint /basis since these just move to conjugate representations. 

Some fundamental work on the Riemann-Hilbert problem was done by Schlesinger [21] . 
He considered the question of constructing new Riemann-Hilbert solutions from a given 
solution RHa(A) = M, in two ways: 

1) Schlesinger examined the fibres of the Riemann-Hilbert map and defined "Schlesinger 
transformations", which move A within the fibres (cf. also JH])- Roughly speaking 
generic fibres are discrete and correspond to certain integer shifts in the eigenvalues of 
the matrices Ai\ geometrically these Schlesinger transformations amount to rational gauge 
transformations with singularities at the poles of the Fuchsian system. 

2) Schlesinger also found how the matrices A can be varied as one moves the pole positions 
a in order to realise the same monodromy data M. (Locally — for small deformations of 
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a — this makes sense as one can use the same loops generating 7ri(P*); globally one should 
drag the loops around with the points a, so on returning a to their initial configuration p 
may have changed by the action of the mapping class group of the m-pointed sphere.) He 
discovered that if the matrices Ai satisfy the following non-linear differential equations, 
now known as the Schlesinger equations, then locally the monodromy data is preserved 
(up to overall conjugation): 



dA, [Ai,A,] _ . ^ dA, ^ [A,,A^] 



(3) ^ = LJ^ if J, and 7^--^ • 

In the generic case such an "isomonodromic deformation" necessarily satisfies these 
equations (up to conjugation). This gives a hint at the difficulty of the Riemann-Hilbert 
problem: even if one knows a solution for some configuration of pole positions, one must 
integrate some nonlinear differential equations to obtain solutions for a deformed config- 
uration. 

This also gives a hint at how one might find some interesting solutions to the Riemann- 
Hilbert problem. Namely since one can move the pole positions one may consider degen- 
erations into systems with fewer poles (for which the problem should be easier). Using 
solutions to these degenerate Riemann-Hilbert problems one can get asymptotics for the 
original solution to the Schlesinger equations and in good circumstances this enables com- 
putation of the solution. This is in effect what we will do below (using the analysis of the 
degenerations in part H and [TK]). 

Suppose we fix an irreducible representation p G Hom(7ri(P*), G). Let Cj C G be the 
conjugacy class containing Mj = p(7i) which we will suppose for simplicity is regular 
semisimple, although this is not strictly necessary. (We are thus considering "generic" 
representations. ) 

Since p is irreducible we know Pj there exists some Riemann-Hilbert solution RHa(A) = 
M. Let C be the adjoint orbit of Ai (in the Lie algebra of n x n complex matrices). 
By genericity we know exp(27rA/— lOj) = Cj. Indeed if in the Riemann-Hilbert map we 
restrict to Ai G Oi then one has Mj G Ci. Also, as mentioned above, the map is equivariant 
under diagonal conjugation and so there is a "reduced Riemann-Hilbert map": 



(4) O := 01 X ■ ■ ■ X Om/IG ^ Ci ® ■ ■ ■ ® CrrJiG =: C 

where the space O is the quotient of {(^i, . . . , Am) | G Oj, ^ = O} by overall conju- 
gation by G and the space C is the quotient of {(Mi, . . . , Mm) \ Mi & Ci^ Mm ■ ■ ■ Mi = l} 
by overall conjugation by G. Generally this map z/a is an injective holomorphic symplectic 
map between complex symplectic manifolds of the same dimension. 

The simplest case is when the representation is rigid, i.e. when the expected dimensions 
of both sides of is zero. Then one knows the RHS of (jH) consists of precisely one point 
and the LHS (at most) one point. 

Our basic strategy is to look at the next simplest case, with the aim of degenerating into 
the rigid case. Since the spaces are symplectic, this corresponds to complex dimension 
two, i.e. both sides of (jU are complex surfaces. 

The principal example of such "minimally non-rigid" systems occurs if we look at rank 
two systems with four poles on the sphere (i.e. n = 2,m = A). Without loss of generality 
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(by tensoring by logarithmic connections on line-bundles) one can work with G = SL2(C) 
rather than GL2(C) and, using automorphisms of the sphere we can fix three of the poles 
at 0,1,00 and label the remaining pole position t. Thus we are considering systems of 
the form: 

By convention we denote the eigenvalues of Ai by ±^j/2 for i = 1,2,3,4. Schlesinger's 
equations imply that the residue A4 = — Ai at infinity remains fixed; we will conjugate 
the system so that = ^ diag(6'4, —O4). The remaining conjugation freedom is then just 
conjugation by the one-dimensional torus T := diag(a, 1/a), a G C*; the space of such 
systems is then three dimensional (quotienting by T yields the surface O) . 

Following pp. 443-446 one may choose certain coordinates x, y, k on this space of 
systems and write down what Schlesinger's equations become. One obtains a pair of 
coupled first-order nonlinear differential equations in x, y (not dependent on k) and an 
equation for k of the form ^ = f(y,t)k. The coordinate k corresponds to the torus 
action, which we can forget about since we are happy to consider Fuchsian systems up to 
conjugation. Eliminating x from the coupled system yields the sixth Painleve equation: 

1 /^l ^ 1 , 1 ^ fdyV fl , 1 , I \ dy 



df^ 2 \y y-1 y-tj\dtj \t t-1 y - 1 J dt 



f-{t-ir V y- (y-ir (y-ty 

where the constants a,/?, 7, 5 are related to the ^-parameters as follows: 

(6) a = (04-1)72, P = -9l/2, 7 = 01/2, S=il-9l)/2. 

Since we will want to go back from a solution of Pyi to an explicit isomonodromic family 
of Fuchsian systems, we will give the explicit formulae for the matrix entries of the system 
in terms of y,y', in appendix 1X1 

Now the bad news is that most solutions to Pyi cannot be written in terms of classical 
special functions. From Watanabe's work j2H| one knows that either a solution is non- 
classical or it is a Riccati solution (corresponding to a reducible or rigid monodromy 
representation p) or the solution y{t) is an algebraic function. 

Since we are interested in explicit solutions corresponding to irreducible non-rigid rep- 
resentations, the only possibility is to seek algebraic solutions to Pyi, in other words 
solutions defined implicitly by equations of the form 

F{y,t) = 

for polynomials F in two variables. We can rephrase this more geometrically: 

Definition 2. An algebraic solution of Pyi consists of a triple {Il,y,t) where 11 is a 
compact (possibly singular) algebraic curve and y, t are rational functions on 11 such that: 

• t : n is a Belyi map (i.e. t expresses 11 as a branched cover of which only 
ramifies over 0, 1, 00), and 

• Using t as a local coordinate on 11 away from ramification points, y{t) should solve 
Pvi, for some value of the parameters a, /5, 7, 5. 
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Indeed given an algebraic solution in the form F{y, t) = one may take 11 to be the 
closure in of the affine plane curve defined by F . That t is a Belyi map on 11 follows 
from the Painleve property of Pyi: solutions will only branch at t = 0, 1, oo and all other 
singularities are just poles. The reason we prefer this reformulation is that often the 
polynomial F is quite complicated and parameterisations of the plane curve defined by F 
are usually simpler to write down. (The polynomial F can be recovered as the minimal 
polynomial of y over C(t), since C{y,t) is a finite extension of C(t).) 

We will say the solution curve 11 is 'minimal' or an 'efficient parameterisation' if y 
generates the field of rational functions on 11, over C(t), so that y and t are not pulled 
back from another curve covered by 11 (i.e. that 11 is birational to the curve defined by 

The main invariants of an algebraic solution are the genus of the (minimal) Painleve 
curve n and the degree of the corresponding Belyi map t (the number of branches the 
solution has over the t-line). 

Now the basic question is: what representations p can we start with in order to obtain an 
algebraic solution to Pyi? Well, the solution must have only a finite number of branches 
and so we can start by looking for finite branching solutions, and hope to prove in each 
case that the solution is actually algebraic. 

The important point is that one can read off the branching of the solution y a.st moves 
around loops in the three-punctured sphere \ {0, 1, oo} in terms of the corresponding 
linear representations p. One finds (cf. e.g. ^ section 4) that p transforms according 
to the natural action of the pure mapping class group (which is isomorphic to vri(P^ \ 
{0, 1, oo}) and thus to the free group on two-letters JF2). Explicitly the generators wi, W2 of 
J^2 act on the monodromy matrices M via Wi = uf where uji fixes Mj for j ^ + (1 ^ 
i ^ 4) and 

(7) uji{Mi, M,+i) = (M,+i, M,+iMiM^\). 

(Incidentally the geometric origins of this in the context of Pyi can be traced back at least 
to Malgrange's work on the global properties of the Schlesinger equations.) The full 
classification of the representations p living in finite orbits of this action is still open, but 
there are some obvious ones: namely if p takes values in a finite subgroup of SL2(C) then 
the JF2 orbit will clearly be finite. 

Thus the program is to take such a finite subgroup F C G, and go through the possible 
representations p : 7ri(P*) V (whose image generates P say) and find the corresponding 
Pvi solutions. The two main problems to overcome in completing this program are: 

1) There are lots of such representations (even up to conjugation), for example for the 
binary tetrahedral group from one knows there are 520 conjugacy classes of triples of 
generators. 

2) We still need to find the Pyi solution explicitly. 

For 1) we proceed as in [5]; by using Okamoto's affine F4 symmetry group of Pyi we can 
drastically reduce the number of classes that arise for each group. It is worth emphasising 
that upon applying an Okamoto transformation the monodromy group may well become 
infinite, and currently there are very few examples of algebraic solutions to Pyi which 
are not equivalent to (or simple deformations of) a solution with finite linear monodromy 
group (see the final remark of sectional below). 
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For 2) we use Jimbo's asymptotic formula (see ^S] and the corrected version in [7j 
Theorem 4). By looking at the degeneration of the Fuchsian system into systems with 
only three poles (hypergeometric systems) and using explicit solutions of their Riemann- 
Hilbert problems, Jimbo found explicit formulae for the leading term in the asymptotic 
expansion of Pyi solutions at zero. Using the Pyi equation these leading terms determine 
the Puiseux expansions of each branch of the solution at zero and, taking sufficiently 
many terms, these enable us to find the solution completely if it is algebraic. 

Philosophically the author views this work as an illustration of the utility of Jimbo's 
asymptotic formula. An alternative method of constructing solutions of Painleve VI has 
been proposed by Kitaev (and Andreev) [201 HI CH] who call it the "RS" method (see 
also Doran for similar ideas, section El below and also for closely related ideas 
of F. Klein). Kitaev [20] conjectures that all algebraic solutions arise in this way and, 
with Andreev, has found some solutions essentially by starting to enumerate all suitable 
rational maps along which a hypergeometric system may be pulled back. 

One of our original aims was to try to ascertain what algebraic Pyi solutions are known, 
up to equivalence under Okamoto transformations and simple deformation (cf. e.g. [Hj 
Remark 15). In other words the aim was to see how much is known of what might be 
called the "non-abelian Schwarz list", viewing Pyi as the simplest non-abelian Gauss- 
Manin connection. The result is that, so far, all the algebraic solutions the author has 
seen have turned out to be related to a finite subgroup of SL2(C) or to the 237 triangle 
group (see section |21 below). ^ As an illustrative example of what can happen consider 
solution 4.1.7.B of jj: At first glance we see t is a degree 8 function of the parameter s 
and so one imagines a solution with 8 branches (and wonders if it is related to one of the 
eight-branch solutions of [S] or ^2] or of section 01 below) . However one easily confirms 
that in fact 



y = y2i = t + 



so it really only has two branches (it was inefficiently parameterised) . In turn one finds 
(for any value of the constant r]oo) this is equivalent to the well-known solution y = \ft. 

On the other hand Jimbo's formula gives us great control, in that we can often go 
directly from a linear representation p to the corresponding Pyi solution. In particular 
the mapping class group orbit of p tells us a priori the number of branches (and lots more) 
that the solution will have. At some point the author realised (see the introduction to [5J) 
that there should be more solutions related to the symmetries of the Platonic solids than 
had already appeared; we have found it to be more efficient to first ascertain directly what 
solutions arise in this way, than for example to enumerate rational maps. (The author's 
understanding is that a theorem of Klein implies that the solutions of sections [HI and HI 
below and of [S] will arise via rational pullbacks of a hypergeometric system, but it is 
not clear if the enumeration started in JTj would ever have found all the corresponding 
rational maps independently.) 



One might be so bold as to conjecture that there are no others, simply because no others have yet 
been seen, in spite of the variety of approaches used. 
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3. The tetrahedral solutions 

In this section we will classify the solutions to Pyi having linear monodromy group 
equal to the binary tetrahedral group T C G = SL2(C). The procedure is similar to that 
used in j2] for the icosahedral group. 

First we examine (as in [5^ section 2) the set S of G-conjugacy classes of triples of gener- 
ators (Ml, M2, M3) of r (i.e. two triples are identified if they are related by conjugating by 
an element of G). (Equivalently this is the set of conjugacy classes of representations p of 
the fundamental group of the four-punctured sphere into F, once we choose a suitable set 
of generators.) From Hall's formulae pi] one knows there are 12480 triples of generators 
of F and dividing by 24 (the size of the image in PSL2(C) of the normaliser of F in G) we 
find that S has cardinality 520. Then we quotient S further by the relation of geometric 
equivalence (cf. ^.5J section 4): two representations are identified if they are related by 
the full mapping class group, or by the set of even sign changes of the four monodromy 
matrices Mj (with M4 = (M3M2Mi)^^). One finds there are precisely six such geometric 
equivalence classes, and by Lemma 9 of 5j this implies there are at most six solutions 
to Pvi with tetrahedral monodromy which are inequivalent under Okamoto's affine F4 
action. 

On the other hand we can look at the set of 6'-parameters corresponding to the represen- 
tations in S. Since Okamoto transformations act by the standard WaiF^) action on the 
space of parameters, it is easy to find the set of inequivalent parameters that arise from 
S, cf. [5^ section 3. (Since they are real we can map them all into the closure of a chosen 
alcove.) We find there are exactly six sets of inequivalent parameters that arise and so 
there are at least six inequivalent tetrahedral solutions. Combining with the previous 
paragraph we thus see there are precisely six inequivalent tetrahedral solutions to Pyi. 

Various data about the six classes and the corresponding Pyi solutions are listed in 
tables 1 and 2. Table 2 lists a representative set of ^-parameters for each class together 
with numbers aij which uniquely determine a triple Mi, M2, M3 in S (and thus the linear 
representation p) for that class with the given 6 values, via the formula 

Tr{MiMj) =2cos{7caij). 





Degree 


Genus 


Walls 


Type 


Alcove Point 


n 


Nonlinear Group 


Partitions 


1 


1 





2 




35, 15, 15, 5 


96 


1 




2 


1 





3 




30, 10, 10, 10 


32 


1 




3 


2 





3 




50, 10, 10, 10 


48 


^2 


1,2 


4 


3 





3 


b^+ 


40, 0, 0, 


72 


53 


3,2 


5 


4 





2 


ab^ 


45, 5, 5, 5 


128 




3 


6 


6 





3 




50, 10, 0, 


144 


A, 


2^3^ 



Table 1 . Properties of the tetrahedral solutions. 



The first two columns of table 1 list the degree and genus of the Pyi solution. The column 
labelled "Walls" lists the number of affine F4 reflection hyperplanes the parameters of the 
solution lie on. The type of the solution enables us to see at a glance which class a given 
element of S lies in: Given Mi,M2,M3,M4 G F their images in PSL2(C) = S03(C) are 
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(n o o n \ 

"2; "3) "4j 


((Ti2, C723, (Jis) 


1 


1/2, 0, 1/3, 1/3 


1/2, 1/3, 1/3 


2 


1/3, 0, 1/3, 1/3 


1/3, 1/3, 1/3 


3 


1/3, 2/3, 2/3, 2/3 


1/2, 1/3, 1/2 


4 


2/3, 1/3, 1/3, 2/3 


1/2, 1/3, 1/2 


5 


1/3, 1/3, 1/3, 1/2 


1/3, 2/3, 1/3 


6 


1/2, 1/3, 1/3, 1/2 


1/3, 1/2, 1/3 



Table 2. Representative parameters for the tetrahedral solutions 



real rotations and we write an "a" for each rotation by half of a turn, a "6" for each 
rotation by a third of a turn, and write nothing for each trivial rotation thus obtained. 
This distinguishes all classes except 3 and 4 which both correspond to four rotations by 
a third of a turn: each Mj thus has parameter = 1/3 or 6*4 = 2/3. For class 3 there are 
always an odd number of each type of (1/3 or 2/3) so we write a minus, and for class 4 
there are always an even number of each type, so we write a plus. 

Finally the rest of table 1 lists the corresponding alcove point (scaled by 60), the number 
n of elements of 5* belonging to each class, the monodromy group of the cover t : 11 — > 
and the unordered collection of sets of ramification indices of this cover over t = 0, 1, 00 
(repeating the last set of indices until three are obtained). Thus for example each solution 
corresponding to row 6 has indices (3,3) over two points amongst {0,1, cxd} and indices 
(1, 1,2,2) over the third. 

All of the tetrahedral solutions have genus zero so we may take 11 to be with param- 
eter s and write the solutions as functions of s. As in the icosahedral case the solutions 
with at most 4 branches are closely related to known solutions. For classes 1 and 2 one of 
the monodromy matrices is projectively trivial and so these rows correspond to pairs of 
generators of the tetrahedral group, i.e. to the two tetrahedral entries on Schwarz's list of 
algebraic hypergeometric functions. The corresponding Pyi solutions are both just y = t 
with the parameters as listed in table 2. As in [5^ one finds class 3 contains the solution 
y = ±-\/t (with the parameters as listed in table 2). Class 4 contains the tetrahedral 
solution 

^ (^^-l)(^ + 2) ^ (^-1)^(^ + 2) 

^ s{s + l) ' (s + l)2(,_2) 

on p. 592 of fT^ (with the parameters as listed in table 2) and is equivalent to a solution 
found independently by Dubrovin [0] (E.31). Also class 5 contains a simple deformation 
of the four-branch dihedral solution in section 6.1 of [12]: 

y-s^ + s + V ^~ 2s + l ' 

that is, this solution is tetrahedral if we use the parameters in table 2, rather than the 
parameters 1/2, 1/2, 1/2, 1/2 for which it is dihedral. 

Thus we are left with one solution, corresponding to row 6. Using Jimbo's asymptotic 
formula to compute the Puiseux expansions etc. (as in [7] section 5, especially p. 193) we 
find the following solution in this class: 
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Tetrahedral solution 6, 6 branches (^^i, 6*2, 6*3, 6*4) = (1/2, 1/3, 1/3, 1/2): 

3(, + 3)(,_l)2' _!)(, + 3)^ 

(We have recently learnt that this is equivalent to solution 4.1.1A in |lj.) It is now easy 
to write down the explicit isomonodromic family of Fuchsian systems in this case, thereby 
solving the Riemann-Hilbert problem for this class of representations p, for an arbitrary 
configuration of the four pole positions (up to automorphisms of P'^). (We will leave for 
the reader the analogous substitutions for the other solutions below.) Using the formulae 
in appendix E] one finds the family of systems parameterised by s G is (up to overall 
conjugation): 



d f Ai Ao A-, 

+ 7T^ + 



where 



dz \ z z — t{s) z — 1 

^ f (s2+3)(^s8-51s"'+99s2-8l) 4s(s-*-9) \ / ^ 

~ I 4(5s6-75s'*+135s2-8l)s(s4_9) -(^s'^j^^'^^^s'^-^Is^+'^Qs^-^i^ j / 

_ f 4(s+3)(s-l)^s2(^s^-s2+3s+9) -2(s+'i)(s-lf[s^+2s+'i) \ I ^ 

~ I -2(s+3)(s-l)^(s3-3s2_9^_9^^5^5_5g4_45^_27^ -4:s'^{s+?.){s-lf[s-^-s'^+?.s+9) j / 

A3 = diag(-^4, 0i)/2 -A1-A2, A = -36 (s^ + 3) (s^ - l)^ [s^ - 9) . 
Note that if the denominator A is zero then t G {0, 1, 00} since 

^ ^_o (^^ + 3)^(^^-3) 
{s + ?,f{s-lf 

Thus the system is well defined for all s in i(:^^(P^ \ {0, 1, 00}) except possibly at s = 00 
(where t = —1). However writing s = 1/s' it is easy to conjugate the system to one well- 
defined also at s = 00. Thus one never encounters configurations requiring a nontrivial 
bundle; the Malgrange divisor is trivial in this situation (in spite of the fact the solution 
y does have a pole at s = 00); indeed one knows the corresponding r function (whose 
zeros lying over P^ \ {0, 1, 00} correspond to nontrivial bundles) satisfies: 

, X ^/ fAi A3 \\ , s^ + 6s^ + 3s^-8s=^-9s2-54s-27 , 
which is nonsingular for all s G t^^(P^ \ {0, 1, 00}). 

Remark 3. Sometimes one is interested in Fuchsian equations with given monodromy, 
rather than systems. To obtain these one may choose a cyclic vector, or more simply 
substitute the Pyi solution into the standard formulae for the isomonodromic family of 
Fuchsian equations. In the present case one obtains the equation: 



rf2 d 

dz^ dz 



where Oi, 02 are respectively: 



12 2 1 

2I ^ 3(2-1) ^ 3{z-t{s)) ~ z-y{s)' 
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1 /7s6-6s^ + 3s^ + 4s3-63s2-54s-27 {s^-2s + 3)s^ 1\ 

z{z-l) V 18(s + 3)3(s- l)3(z-t(s)) 9{s + 3){s - iy{z - y{s)) ~ 18 J ' 

For generic s this is a Fuchsian equation with non-apparent singularities at z = 0, l,t, oo 
and an apparent singularity ai z = y, realising the given (projective) monodromy repre- 
sentation. In special cases (when y = O,l,t,oo) it will have just the four non-apparent 
singularities (and will thus be a so-called "Heun equation"). For example specialising to 
s = one finds y = 0,t = —1 and the equation becomes that with 

_ 1 2 2 _ 1 

~ ~2i^ 3(;2-l) ^ 3(^+1) ~ ~ 18(^- 1)(2 + 1) 

which is a Heun equation whose projective monodromy representation is that specified 
by row 6 of table 2. 

Remark 4. At the editor's request we will explain how one may verify directly that these 
Pvi solutions actually do correspond to Fuchsian systems with linear monodromy repre- 
sentations as specified by table 2. For the rigid cases, rows 1 and 2, this is immediate, by 
rigidity. For the others, first one may check that the solutions actually do solve Pyi- This 
can be done directly (by computing the derivatives of y with respect to t and substitut- 
ing into the Pyi equation).^ Having done this we know the formulae of appendix \K\ do 
indeed give an isomonodromic family of Fuchsian systems. To see it has the monodromy 
representation specified by table 2 we first compute the Puiseux expansions at of each 
branch of the function y{t) (only the leading terms will be needed). On the other hand 
Jimbo's asymptotic formula (in the form in [Zj Theorem 4) computes the leading term 
in the asymptotic expansion of the Pyi solution corresponding to the given monodromy 
representation p (the leading term is of the form at'^ where a and b are explicit functions 
of 6*1, 6^2, 6^3, 6*4, (Ji2, (723, (Jis). Then it is sufficient to check this leading term equals one 
of the leading terms of the Puiseux expansions of y{t). The logic is that, in the cases 
at hand, the leading term determines the whole Puiseux expansion (using the recursion 
determined by the Pyi equation) and this is convergent so determines the solution locally, 
and thus globally by analytic continuation. (For the solutions we construct here this is 
automatic since we constructed the solution starting with the results of Jimbo's formula.) 
Some simpler, but not entirely conclusive, checks follows: 

1) compare the monodromy of the Belyi map t with the JF2 action ((Zj) on the conjugacy 
class of the representation p (if we didn't know better it would appear as a miracle that 
the solution, constructed out of just the Puiseux expansion at 0, turns out to have the 
right branching at 1 and 00 too). 

2) compute directly the Galois group of one of the Fuchsian systems in the isomon- 
odromic family. (Together with the exponents this goes a long way to pinning down the 
monodromy representation.) There are various ways to do this, one of which is to convert 
the system into an equation (e.g. via a cyclic vector) and use the facility on Manuel 
Bronstein's webpage: 

http://www-sop.inria.fr/cafe/Manuel.Bronstein/sumit/bernina_demo.html 

^ To aid the reader interested in examining the solutions of this article (and to help avoid typographical 
errors) a Maple text file of the solutions has been included with the source file of the preprint version 
on the math arxiv (|math.DG/0501464||. This may be downloaded by clicking on "Other formats" and 
unpacked with the commands 'gunzip 0501464. tar' and 'tar -xvf 0501464. tar', at least on a Unix system. 
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(This requires finding a suitable rational point on the Painleve curve, which, if possible, 
is easy in the genus zero cases, and not too difficult using Magma in the genus one cases.) 

4. The octahedral solutions 

For the octahedral group we do better and find more new solutions. In this case, by ^T] 
or direct computation, S has size 3360, which reduces to just thirteen classes under either 
geometric or parameter equivalence. Thus there are exactly thirteen octahedral solutions 
to Pvi, up to equivalence under Okamoto's affine F4 action. 

Data about these classes are listed in tables 3 and 4. In this case the type of the solution 
may contain the symbol "(7" which indicates that one of the corresponding rotations in SO3 
is a rotation by a quarter of a turn. Also, in some cases rather than list the monodromy 
group of the cover t : 11 — we just give its size. 





Degree 


Genus 


Walls 


Type 


Alcove Point 


n 


Group (size) 


Partitions 


1 


1 





1 


abg 


(65, 35, 25, 5)/2 


192 


1 




2 


1 





2 




25, 10, 10, 5 


96 


1 




3 


2 





2 




45, 15, 10, 10 


96 




1,2 


4 


3 





1 


abg'^ 


40, 10, 5, 5 


288 


^3 


3,2 


5 


4 





2 


ag^ 


(75, 15, 15, 15)/2 


128 


A, 


3 


6 


4 





3 


g' 


30, 0, 0, 


32 


A, 


3 


7 


6 





2 


d^hg 


(95, 25, 5, 5)/2 


576 


24 


2^3^2 4 


8 


6 





2 




35, 5, 0, 


288 


36 


3,2 4 


9 


8 





1 


ab'^g 


(85, 15, 15, 5)/2 


768 


576 


2^3,2M 


10 


8 





3 


a'g' 


45, 15, 0, 


192 


192 


3^2 3^ 


11 


12 





3 




50, 10, 0, 


288 


576 


2^3^2M^ 


12 


12 


1 


3 






288 


96 




13 


16 





3 


a-^g 


(105, 15, 15, 15)/2 


128 


3072 


2^34 



Table 3. Properties of the octahedral solutions. 



The octahedral solutions with at most 4 branches correspond to the following known 
solutions. As in ,5] one finds: The first two classes correspond to the octahedral entries 
on Schwarz's list of algebraic hypergeometric functions (and the Pyi solution isy = t with 
the parameters indicated in table 4). Solution 3 is y = ±\/t with the parameters listed in 
table 4, solution 4 has 3 branches and is a simple deformation of the 3-branch tetrahedral 
solution above (namely it is the solution in equation (jH)), but with the parameters given 
in table 4), solution 5 is a simple deformation of the 4-branch dihedral solution (namely 
it is the solution in equation Q, but with the parameters given in table 4), and solution 
6 is the 4-branch octahedral solution 

y = 7(^y ' = s^is - 2) 

on p. 588 of [13 , with the parameters as in table 4, which is equivalent to a solution found 
independently by Dubrovin P (E.29). 

For the remaining 7 solutions, rows 7-13, we will construct an explicit solution in each 
class using Jimbo's asymptotic formula. More computational details appear in appendix 
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/'a ZD ZD N 
[Ui, ^2, t^s, t:74j 


((712, 0-23, CTis) 


1 
i 


I/O n 1 /Q 1 //I 
i//:, U, i/6, i/4 


I/O 1 /Q 1 //I 

i/2, i/3, i/4 


o 
/ 


1 /Q n 1 //I 1 //I 
i/d, U, 1/4, 1/4 


1 /Q 1 //I 1 //I 

1/d, 1/4, 1/4 


q 


1 /Q 1 //I 1 //I O /Q 

1/3, 1/4, 1/4, 2/3 


I/O I/O I/O 

i/2, i/2, 1/2 


A 

4 


I/O 1 //I 1 //I 9 /Q 

J-/^, -1-/4, J-/4, Z/6 


1/0 1 /Q Q //I 

1/2, 1/d, d/4 





1 //I 1 //I 1 //I I/O 

-1-/4, i/4, 1/4, 1//: 


1 /Q I/O 1 /Q 

1/3, 1/2, 1/3 


D 


1 //I 1 //I 1 //I 1 //I 
-1-/4, i/4, i/4, i/4 


1 /Q n 1 /Q 
i/d, U, i/d 


7 


1/9 1/9 1 /zL 9 /'^ 

-L/-^! -L/^) -L/^j 2/<J 


1/9 1/9 1 /'^ 
-1-/2, 1/2, 1/J 


8 


1/3, 3/4, 1/3, 3/4 


1/2, 3/4, 1/3 


9 


1/3, 1/4, 1/2, 2/3 


1/2, 2/3, 3/4 


10 


1/2, 1/4, 1/2, 3/4 


2/3, 2/3, 1 


11 


1/3, 1/2, 1/2, 2/3 


1/2, 1/2, 1/4 


12 


1/2, 1/2, 1/2, 2/3 


1/2, 1/4, 2/3 


13 


1/2, 1/2, 1/2, 3/4 


1/2, 2/3, 1/3 



Table 4. Representative parameters for the octahedral solutions 



O (We have recently learnt that solutions 8 and 10 are equivalent to those of ^H], 3.3.3 
top of p. 22, and 3.3.5 bottom of p. 23, respectively.) The formulae obtained are as follows: 

Octahedral solution 7, 6 branches (^i, 62, 63, 6^) = (1/2, 1/2, 1/4, 2/3): 
9s(2s3-3s + 4) _ 27^2 

^~ 4(s + l)(s- l)'(2s2 + 6s + l)' ~ 4{s^-lf 

Octahedral solution 8, 6 branches (^i, 62, 63, 64) = (1/3, 3/4, 1/3, 3/4): 
(2s2-l)(3s-l) (3s -1)' 



y 



2s (2 s2 + 2 s - 1) (s - 1) ' 8 (2 s2 + 2 s - 1) (s - 1) 



Octahedral solution 9, 8 branches (^i, ^2, ^3, ^4) = (1/3, 1/4, 1/2,2/3): 
_ s3(2s2-4s + 3)(s2-2s + 2) _ fs^2s^ -4s + 3) 



(2s2-2s + l)(3s2-4s + 2) ' V 3s2-4s + 2 

Octahedral solution 10, 8 branches (^1,^2,^3,^4) = (1/2,1/4,1/2,3/4): 
32s(s + l)(5s2 + 6s- 3) 10245^ (s + 1)^ 



(s2 + 2 s + 5) (3 s2 + 2 s + 3)^ ' (s^ + 6 s + 1) (3 + 2 s + 3)^ 

Octahedral solution 11, 12 branches (^i, ^2, ^3, ^4) = (1/3, 1/2, 1/2,2/3): 
_ (s + 1) (7 + 16 + 4 s2 - 4) r ^ _ f (s + if r 



^ s3(s-2)(s4-4s2 + 32s-28) ' V " 2) s 

where r = 4 (3 - 4 s + 2) / (s^ + 4 s + 6). 
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The next solution, number 12, has genus one. In this case we take 11 to be the elhptic 
curve 

= (2 s + 1) (9 + 2 s + 1) . 
As functions on this curve the solution is: 

Octahedral solution 12, genus one, 12 branches (^^i, 6*2, 6*3, 6*4) = (1/2,1/2,1/2,2/3): 

_ 1 45 + 20 + 95 + 92 + 39 - 3 

^~ 2 ^ 4 (5s2 + l) {s + lfu 

_ 1 s (2 s + 1)^ (27 + 28 + 26 s2 + 12 s + 3) 

~ 2 ^ {s + lfu^ 

Finally the last solution, number 13, has 16 branches and genus zero. This is possibly 
the highest degree genus zero solution amongst all algebraic solutions of Pyi. It is also 
special since it has no real branches. Thus necessarily the parameterisation of the solution 
is not defined over Q although the solution curve F{y, t) = itself has Q coefficients, as 
does the Belyi map t. 

Octahedral solution 13, 16 branches, (^1,^2,^3,^4) = (1/2, 1/2, 1/2,3/4): 

_ (1 + z) {s^ - 1) {s^ + 2ts+l){s^ -2is + if c _ {s^-lf{s^ + 6s^ + lf 
^~ 4 s (s2 + i) (s2 - if (s2 + (1 + i)s -i)d ' ~ 32 s^s^ + if 
where 

c = s^-(2- 2i) - (6 + 2i) s^+{10 + 2i) + A is^ + {10- 2i) + (6 - 2i) - (2 + 2i) s - 1 

= / - (3 + 3i) + 3 + (4 - Ai) + 3 + (3 + 3i) s + z 

Remark 5. The author has recently understood that an alternative way to construct some 
(but not all) of these tetrahedral and octahedral solutions would have been to use the 
quadratic transformations of ^H]- For example tetrahedral solution 6 could have been 
obtained from tetrahedral solution 4 in this way (a fact that was apparently not noticed 
in PJ). It is debatable whether this would have been simpler for us than the direct method 
used here, given what had already been done in [71IS]. (The quadratic transformations 
were crucial however to construct the higher genus icosahedral solutions, cf. [^.) 

5. Infinite monodromy groups 

In this final section we will give some examples of solutions corresponding to non-rigid 
representations p into some infinite subgroups of G = SL2(C). The point is that the 
method we are using to construct Pyi solutions should work provided that the JF2 orbit 
of p is finite, and above we just used the finiteness of the image of p as a convenient way 
of ensuring this. 

Thus we are looking for representations p having finite JF2 orbits, or more concretely, 
matrices Mi,M2,M3 G G having finite orbit under the action ((7j). (Such an JF2 orbit, 
on conjugacy classes of representations, gives the permutation representation of the Belyi 
cover t : n — for the corresponding Pyi solution. We would like to find interesting 
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Ti orbits in order to find interesting Pyi solutions.) So far tliere appear to be four ways 
(apart from guessing) of finding representations p liaving finite T2 orbits. 

Firstly one can just set the parameters to be sufficiently irrational in one of the families 
of solutions. (For example y = \/t is a solution provided 9i + 64 = 1,62 = O3 amongst 
other possibilities.) 

Secondly one can sometimes apply an Okamoto transformation to a known solution and 
change p into a representation having infinite image. For example if we take the 16 branch 
octahedral solution above and apply the Okamoto transformation corresponding to the 
central node of the extended D4 Dynkin diagram, then we obtain a Pyi solution whose 
corresponding linear representation has image equal to the (2, 3, 8) triangle group. To 
see this we recall [TU [7j that Okamoto's affine D4 action does not change the quadratic 
functions Tr(MiMj) = 2cos{naij) of the monodromy data, only the ^-parameters. In this 
case the 16 branch octahedral solution has data 

6 = (1/2, 1/2, 1/2, 3/4), a = (1/2, 2/3, 1/3) 

on one branch and the solution after applying the transformation has data 

9 = (3/8, 3/8, 3/8, 5/8), a = (1/2, 2/3, 1/3). 

One may show that the image in PSL2(C) of the corresponding triple (Mi, M2, M3) gen- 
erate a (2,3,8) triangle group (see appendix lB|l . The corresponding solution to Pyi is 
given by the formula 

Z/238(s) = y{s) + — — 

where y, t, 9i are as for the 16 branch octahedral solution and p as in appendix ^ Ex- 
plicitly one finds the solution is 

2,3,8 solution, genus zero, 16 branches, (6'i, 6^2, 6'3, 6*4) = (3/8,3/8,3/8,5/8): 

_ (1 + i) {s^ - 1) {s^ + 2is + l){s^ -2is + ifd' 



with t and d{s) as for the 16 branch octahedral solution, and d'{s) = d{s). In turn, via 
the formulae in appendix this yields the explicit family of Fuchsian systems having 
projective monodromy group the (2, 3, 8) triangle group. 

Thirdly, the idea of [Zj was to use a different realisation of Pyi as controlling isomon- 
odromic deformations of certain 3x3 systems. The corresponding monodromy represen- 
tations were subgroups of GL3(C) generated by complex reflections, and again one will 
obtain finite branching solutions by taking representations into a finite group generated 
by complex reflections. Applying this to the Klein complex reflection group led to an 
algebraic solution to Pyi with 7 branches. Moreover [7j described explicitly how to go 
between this 3x3 picture and the standard 2x2 picture used here, both on the level 
of systems and monodromy data. The upshot is that if we substitute the Klein solution 
into the formulae of appendix El below, then we obtain an isomonodromic family of 2 x 2 
Fuchsian systems with monodromy data on one branch given by: 

(11) 9 = (2/7, 2/7, 2/7, 4/7), a = (1/2, 1/3, 1/2). 

This determines a representation into (a lift to G of) the (2, 3, 7) triangle group (and one 
may show as in appendix El its image is not a proper subgroup). Moreover it was proved 
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in [7j that this cannot be obtained by Okamoto transformations from a representation 
into a finite subgroup of SL2(C). 

Fourthly one may obtain such representations by pulhng back certain hypergeometric 
systems along certain rational maps (cf. Doran 8J and Kitaev JH])- This is closely related 
Klein's theorem that all second order Fuchsian equations with finite monodromy group 
are pull-backs of hypergeometric equations with finite monodromy. The basic idea is as 
follows. 

Label two copies of by u and d (for upstairs and downstairs). Choose four integers 
no, ^1, Hoc, N > 2 and suppose we have an algebraic family of branched covers 

of degree N, parameterised by a curve 11 say, such that: 

1) TT only branches at four points 0, 1, oo and at a variable point x G P^, 

2) All but four of the ramification indices over 0, 1, oo divide no, ni, Uoo respectively. In 
other words if {cij} are the ramification indices over i = 0, 1, oo then as j varies precisely 
four of the numbers 

no ' ni ' noo 

are not integers. Let t be the cross-ratio of the corresponding four ramification points of 
P^, in some order, and so we have a coordinate on P^ such that these four points occur 
at 0, t, 1, oo. 

3) TT has minimal ramification over x, i.e. tt ramifies at just one point over x, with 
ramification index 2. 

The idea of jH], JBi is to take a hypergeometric system on P^ with projective monodromy 
around i equal to an n^'th root of the identity, for z = 0, 1, oo and pull it back along tt. 
One then obtains an isomonodromic family on P^ with non-apparent singularities at 
0, t, 1, oo and possibly some apparent singularities at the other ramification points. All of 
the apparent singularities can be removed, for example by applying suitable Schlesinger 
transformations, yielding an isomonodromic family of systems of the desired form. 

In particular the problem of constructing algebraic solutions of Pyi now largely becomes 
a purely algebraic problem about families of covers, although it is only conjectural that 
all algebraic solutions arise in this way. 

However the algebraic construction of such covers seems difficult. First one has the 
topological problem of finding such covers, then one needs to find the full family of covers 
explicitly (this amounts to finding a parameterised solution of a large system of algebraic 
equations, typically with one less equation than the number of variables, so the solution 
is a curve). See |^ for some interesting examples however (but one should be aware that 
some of these solutions are equivalent to each other and to known solutions via Okamoto 
transformations) . 

Our perspective here is that just the topology of the cover is enough to determine the 
monodromy of the Fuchsian equation on P^ and we can then apply our previous method 
[7j to construct the explicit Pyi solution. In other words just one topological cover vr gives 
us the desired representation p living in a finite J^2 orbit. 

To find some interesting topological covers we consider the list appearing in Corollary 
4.6 of Here Doran classified the possible ramification indices of the cover vr in the cases 
where the monodromy group of the hypergeometric system downstairs is an arithmetic 
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triangle group in SL2(M). (Contrary to the wording in [5] this does not determine the 
topology of the cover.) We will content ourselves with looking at the last four entries of 
Doran's list, which say that the integers N,no,ni,noo and the ramification indices are: 



10, 


2,3,7 


[2,.. 


.,2], [3,3,3,1], [7,1,1,1] 


12, 


2,3,7 


[2,.. 


.,2], [3, 3, 3, 3], [7, 2, 1,1,1] 


12, 


2,3,8 


[2,.. 


.,2], [3,3,3,3], [8,1,1,1,1] 


18, 


2,3,7 


[2,.. 


.,2], [3,. ..,3], [7,7,1,1, 1,1] 



The basic problem now is to find such covers topologically, in other words to find 
the possible permutation representations. (The cover of the four punctured sphere \ 
{x, 0,1,00} is determined by its monodromy, which amounts to four elements of Sym^ 
having product equal to the identity and whose conjugacy classes — i.e. cycle types — are 
as specified by the given ramification indices.) 

The simplest way to do this is to draw a picture. Suppose we fix x = —1 and cut 
along the interval / := [—1, 00] from —1 along the positive real axis. Then the preimage 
of / under vr will be a graph in P^ with vertices at each point of 7r~^{{x, 0, 1, 00}). The 
complement of the graph will be the union of exactly connected components which 
are each mapped isomorphically by vr onto P^ \ /, and in particular the boundary of each 
component is the same as the boundary of P^\/. These connected components correspond 
to the branches of vr and the graph specifies how to glue them together. In particular the 
graph determines the permutation representation of the cover, since it shows us how to 
lift loops in the base to paths in P^; we just cross the corresponding edges upstairs, and 
note which connected component we end up in. 

Thus we need to draw the graphs in P^. There are four types of vertices, depending 
on if they lie over —1,0,1, 00, which we could draw as circles, squares, blobs and stars 
(say) respectively. The number of each type of vertices is just the number of points of P^ 
lying over the corresponding point amongst — 1,0, 1,cxd. The corresponding ramification 
indices give the number of edges coming out of each vertex to each of the neighbouring 
vertices, and our task is to join these edges together in a consistent manner. 

For example for the first row of the above list, there should be 10 branches and, by 
examining the ramification indices, we see we need to draw a graph on P^ out of the 
following pieces: 

• 8 circles with 1 edge emanating from each, and 1 circle with 2 edges, 

• 5 squares with 4 edges, 

• 3 blobs with 6 edges and one blob with 2 edges, and 

• 1 star with 7 edges and 3 stars with 1 edge. 

The graph should divide the sphere into 10 pieces and: 

• Each edge from a circle should connect to a square, 

• Two edges from each square should connect to a circle and the other two should 
connect to a blob (and, going around the square, the edges should alternate between 
going to circles and blobs), 

• Similarly half the edges from each blob should connect to squares and, again alter- 
nating, the other half should connect to stars, 

• Finally each edge from a star should connect to one of the blobs. 
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We leave the reader to draw such a graph (there are 15 possibihties).'^ Given any such 
graph we can write down the monodromy of the pulled back Fuchsian system on in 
terms of that of the hypergeometric system downstairs. Here the projective monodromy 
downstairs is a (2, 3, 7) triangle group: 

A237 { a,h,c \ = = c' = cha = I ) 

which can be realised as a subgroup of PSL2(C) in various ways (the standard represen- 
tation into PSL2(M) plus its two Galois conjugates, lying in PSU2). 

Puncture at the four exceptional vertices (namely the 3 stars with 1 edge and the blob 
with 2 edges) and choose generators /i, ■ ■ ■ , /4 of the fundamental group of this punctured 
sphere, with l^o ■ ■ ■ oli contractible. Then we can compute the image under vr of each li 
in P^ \ {0, 1, 00} and thereby write the monodromy of the pulled back system as words in 
a, 6, c G A237. With one such graph we obtained: 

Li = caca^^c^^, L2 = c, L3 = c^^a^^cac, L4 = c~'^bc^ 

where Lj is the projective monodromy around /j. By construction L4 ■ ■ ■ Li = 1 in A237. 
Now by choosing an embedding of A237 in PSL2(C) we get Lj G PSL2(C) and we can 
lift each Li to a matrix Mj G G, (and possibly negate M4 to ensure M4 ■ ■ ■ Mi = 1). We 
obtain the representation p with data 

e = (2/7, 2/7, 2/7, 1/3) a = (1/3, 1/3, 1/7). 

This completes our task of producing a representation in a finite JF2 orbit. Now we can 
apply our previous method to construct the corresponding Pyi solution. Immediately, by 
computing the J-2 orbit of the conjugacy class of p, we find the solution has genus 1 and 
18 branches, and that the parameters are not equivalent to those of any known solution. 
Moreover it turns out that Jimbo's asymptotic formula may be applied to 17 of the 18 
branches, and the asymptotics on the remaining branch may be obtained by the lemma 
in section 7 of jS]. Using this we can get the solution polynomial F explicitly from the 
Puiseux expansions, and then look for a parameterisation of F. The result is: 

2, 3, 7 solution, genus one, 18 branches, (6*1, 62, 63, 64) = (2/7, 2/7, 2/7, 1/3): 

_ 1 (3 - 2 - 4 - 204 - 536 - 1738 - 5064 - 4808 s - 3199) u 
^ ~ 2 4 (36 + 196 + 189 s2 + 756 s + 154) (5^ + 5 + 7) (s + 1) 

_1 (s^ - 84 8*^ - 378 - 1512 - 5208 - 7236 - 8127 s - 784) M 
^~ 2 432s(s + l)'(s2 + s + 7)' 

where 

= si^s"^ + S + 7). 

This solution is noteworthy in that currently there is no known relation to a Fuchsian 
system with finite monodromy group (one might speculate as to the existence of another 
realisation of Pyi in which this solution corresponds to such a Fuchsian system, but this 
is unknown). 

"'To count the possibilities, one may use Theorem 7.2.1 in Serre's book |2S| to count the number of such 
permutation representations, and then divide by conjugation action of the symmetric group, carefully 
computing the stabiliser. To find all possibilities we draw some and then apply the natural action of the 
pure three-string braid group to see if we get them all — here all 15 are braid equivalent. 
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For the other three entries on the excerpt of Doran's hst above, we do not seem to get 
new solutions, but it is interesting to identify them in any case. 

The second entry, a family of degree 12 covers, turns out to give the Klein solution of 
[7j. The explicit family of covers has been found more recently in ^H] P-27. There are 
7 different graphs one could draw, one of which is symmetrical and they are all braid 
equivalent. Using one of these graphs we obtain, as above, the words: 

Choosing an appropriate embedding of A237 and lifting to G we obtain the representation 
p specified in flllj) . In particular this gives a convenient way to prove that the projective 
monodromy group of the family of 2 x 2 Fuchsian systems determined by the Klein solution 
is A237. We just need to show that the Lj generate all of A237, which we will do in appendix 
E] below. 

The third entry indicates a family of degree 12 covers along which one should pull back 
the (2, 3, 8) triangle group. This time there are 7 graphs one could draw but they are not 
all braid equivalent, there are two P3 orbits, distinguished by the fact that the monodromy 
group of the cover is either Sym^g ^ group of order 1536. For the degenerate case one 
finds the Pyi solution has just two branches and is y = t ± — 1) with parameters 
9 = (1, 1, 1, 7)/8. (This is just the transform of the square root solution y = y/i under the 
Okamoto transformation {y,t) (fEf) ttt))- The other case is more interesting; for one 
graph we obtain: 

Li = aca~^, L2 = c~'^ cac^ , L3 = caca~^c~^, L4 = caca'^ a 

where now a, 6, c generate the (2, 3, 8) triangle group: 

A238 = ( a, 6, c I = 1)^ = = cha = l). 

Now we can choose an embedding of A238 into PSL2(C) and a lift to G (negating M4 if 
necessary) such that we obtain the representation p with data 

e = (3/8, 3/8, 3/8, 5/8), a = (1/2, 2/3, 1/3). 

This is precisely that obtained above by applying an Okamoto transformation to the 16 
branch octahedral solution (and gives a convenient way to prove, in appendix^ that the 
projective monodromy group is A238). 

Finally there are 9 graphs corresponding to the last row of Doran's list, all braid equiv- 
alent. Even though the graphs are the most complicated in this case (and there is a quite 
attractive one with 4- fold symmetry), this case leads again to the 2-branch Pyi solution 
y = t ± \/t{t — 1) with the parameters 6 = (1, 1, 1, 6)/7 (and a = (1/2, 1/2, 5/7) on one 
branch) . 

In conclusion we should mention that we do not know any other finite J^2 orbits of 
triples of elements of SL2(C) (e.g. up to isomorphism as 'sets with jF2-action'); so far 
they all either come from a finite subgroup or one of the two (2, 3, 7) cases (the Klein 
solution or the genus one solution above). 

Appendix A. 

Here are the explicit formula from [TUj for the residue matrices Ai, of the isomonodromic 
family of Fuchsian systems corresponding to a Pyi solution y{t) with parameters ^1, . . . , 6^4. 
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The matrix entries are rational functions of y,t,y' = {Oi}- (Our coordinate x is denoted 
z in [Kj and is related simply to p which is the usual dual variable to y = q in the 
Hamiltonian formulation of Pyi- Also one should add diag(6'j, 6'j)/2 to our Ai to obtain 
that of jl6).) 



A, 

where 



E-ki(t + l) , ,E + tBAy -l)xkl-tkik2 

''=y — — ' ^2 = (y-^) ^^31)^^ 

, ..E + 9A{y-t)x-kh-kik2 

= - {F^, ' 

Oi 02 03 ^ 0i + {t- 1) y' , e2-i + y' , e^i-ty' 

X = p , 2p = 1 1 

y y -t y -I y y-t y -l 

y y — t y — 1 



tZi' t{t-l)z2' {t-l)z3 

E = y{y- l){y - t)x^ + {O^iy - t) + t02(y - 1) - 2k2{y - l)(y -t))x + k^y - A:2(03 + ^^2) 
ki = (04 -01-92- 03)/2, k2 = {-04 - 01 - 02 - 03)/2. 



Appendix B. 

Proposition 6. The 2x2 Fuchsian systems corresponding to the Klein solution and to the 
18 branch genus 1 solution of section\^have projective monodromy group isomorphic to 
A237, and that corresponding to the transformation of the 16 branch octahedral appearing 
in section\^has projective monodromy group isomorphic to ^2zs- 

Proof. Since in section |S] the projective monodromy groups were expressed as words 
in the generators of the respective triangle groups, it is sufficient to check in each case 
that these words are in fact generators. To do this we will repeatedly use the fact that in 
the group 

A = ( a, 6, c I a"^ = = = cba = 1 ) 

one has = bc~^ and where in general we write for y ^xy. (These are 

easily verified, for example the first is true since b~^cbc = b^^{cb){cb)b^^ = = b, using 
the fact that a = = cb.) In particular we immediately deduce A = (c, c^) = (c, c^ ). 

Now each case is an easy exercise. For the Klein case we need to show (Lj, i = 1, 2, 3) = 
A where n = 7 and 

-Li — C , 1^2 — C , iv3 — C 

Up to conjugacy in A, we have {L2, L^) = {p, c) where p = c"-^ °- . However, using a = cb 

cb'^ b'^ b-^ J 

we see p = c = c = c so we are done. 
For the other (2, 3, 7) case corresponding to the genus one solution we have 

Li = c'^'="\ L2 = c, Ls = c'''. 

Thus {L2, L3) = (c, c") = (c, c*) since a = cb. 

For the (2, 3, 8) case we have Li = c"", L2 = c""^ , L3 = c"''^ . Up to conjugacy (Li, L3) = 
(c, c"'^ and as in the Klein case above c"'^ " = bcb~^ . □ 
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Appendix C. 

At the request of the editors and of A. Kitaev, we will add some remarks to aid the 
reader interested in reproducing the results of this article. The main results are of two 
types: 1) classification of Pyi solutions coming from the binary tetrahedral and octahedral 
groups and 2) construction of explicit Pyi solutions using Jimbo's asymptotic formula. 
For both 1) and 2) the details are parallel to those described in [5^ and /7j resp., with the 
precise references as in the body of this article. For 1) there are 3 steps: 

• Prove that the relation of Okamoto equivalence is sandwiched between the relations of 
geometric and parametric equivalence, i.e. in symbols one has GE =^ OE =^ PE. The 
second arrow is immediate by definition and the first is proved in Lemma 9 of |5j. 

• Compute the parameter equivalence classes in the set of parameters coming from triples 
of generators of either the tetrahedral or octahedral group. This is as in section 3 of 
[Sj. One first writes down the set of possible parameters 6. This is a finite subset of 

C M^. Then one uses a simple algorithm to move each of these points into a chosen 
affine F4 alcove, using the standard action of the affine Weyl group Wa^F^) on (this 
is entirely standard and the details are written in [Sj Proposition 6). Then we count the 
number of distinct alcove points obtained. By definition this is the number of "parameter 
equivalence classes". 

• Compute the geometric equivalence classes in the set of linear representations p coming 
from either the binary tetrahedral or octahedral group. This amounts to computing the 
orbits of an explicit action of a group on a finite set (of size 520, 3360 resp.) and is 
carefully described in section 4 of 

Some confidence that there is no computational error comes from the fact that the 
geometric and parametric equivalence classes turn out to coincide in both the cases con- 
sidered here. Also Hall's formulae [11^ (computing the number of generating triples) gives 
confidence that all the generating triples have been computed correctly — since we get the 
right number of them. (In principle one can go through all triples of elements of the finite 
group r C SL2(C) and throw out those that do not generate P. In the two cases at hand 
this is feasible, but some simple tricks are useful in the icosahedral case.) 

Now we will move on to 2), constructing the solutions. The main steps of the procedure 
used are as in jjj (see especially p. 193). However with experience various tricks have been 
developed to speed up the computation, so we will also detail some of these below (they 
are inessential if one has a fast enough computer, as presumably future readers will have). 
The underlying strategy is analogous to that used in pi]| although we do not in fact use 
any of their results. (It was particularly troublesome to get the correct form of Jimbo's 
formula in [7j, which is the main ingredient and was not used in J^.) 

The basic steps are as follows: 

1) We start with a linear representation p living in a finite mapping class group orbit. 
The conjugacy class of p is encoded in the seven-tuple 

62, 03, (^A, O'u, 0-23, 0"i3. 

Specifying these seven numbers is equivalent to specifying the numbers rrii = 2cos(7r6'j), 
rriij = 2cos(7r(Tjj) provided we agree Oi,aij G [0, 1]. We compute the orbit of this 7-tuple 
under the pure mapping class group = JF2. The formula for this action is given in |3] 
section 4 (cf. also ((Zj) above). This gives a list, of length N say, of 7-tuples, one for each 
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branch of the corresponding Pyi solution. The values of the 6''s will not vary on different 
branches so the branches are parameterised by the values of the sigmas. Let their values 
on the kth branch be denoted afj, k = 1, . . . ,N. 

2) Plug each 7-tuple into Jimbo's asymptotic formula (in the form in 7] Theorem 4). 
This gives N leading terms yk = cikt^'' + ■ ■ ■ for = 1, . . . , of the Puiseux expansion at 
of the Pvi solution y{t) on the branches. One will have bk = 1 — cr^2 but is given by 
a very complicated, but explicit, formula. (Jimbo's formula is not always applicable — cf. 
the discussion of 'good' solutions in but often there is an equivalent solution for which 
Jimbo's formula can be applied on every branch, or there is a degeneration of Jimbo's 
formula (as in jTUj or j3] Lemma 19) which will compute the remaining leading terms.) 

3) Compute lots of terms in the Puiseux expansions of the solutions y{t) on each branch. 
These will be expansions in s = t^^^'' where dk is the denominator of bk (when written in 
lowest terms) . Geometrically dj. is the number of branches of y that meet the given branch 
over t = 0, i.e. it is the cycle length of the given branch in the permutation representation 
of the solution curve as a cover of the t-line. The expansions are computed recursively by 
substituting back into the Pyi equation; at each step this leads to a linear equation for 
the coefficient of the next term in the expansion. 

4) Use these expansions to determine the coefficients of the solution polynomial F{y, t). 
(This determines y as an algebraic function of t by the condition F{y,t) = 0.) Since F 
is a polynomial (of degree N in y) this is clearly possible since we have arbitrarily many 
terms of each Puiseux expansion; F{yk{s), s'^'') = for all s and for each branch y^ of 
the solution. (Thus in principle just one expansion is needed, not the expansion for all 
branches.) Given F{y,t) one may check symbolically that it specifies a solution to Pyi, 
using implicit differentiation. 

5) Compute a parameterisation of the resulting curve F{y, t) = 0. (As mentioned in the 
acknowledgments the author is grateful to Mark van Hoeij for help with this last step.) 
In general this will be simpler to write down than the polynomial F. 

Now we will list some of the tricks we have found useful in carrying out the above steps. 

1) One needs to convert the numbers ak given by Jimbo's formula into algebraic numbers. 
In the examples so far this can be done by raising (and/or its real/imaginary parts) to 
the dk-th power until a rational number is obtained (which can be ascertained by looking 
at continued fraction expansions). 

2) Reduce the number of Puiseux expansions to compute: the d^ branches which meet 
the given branch over zero will have Galois conjugate expansions. These can be obtained 
from one another by multiplying s by a dk-th root of unity. Also, when choosing which of 
these dk expansions to actually compute it is good to choose the real one, if possible. (Also 
sometimes some expansions are complex conjugate to others so further optimisations are 
possible.) 

3) Reduce the degree of the field extension used to compute the expansion: In computing 
the Puiseux expansions one is often working over a finite extension of Q, such as Q(6^/^). 
Often the degree of this extension can be reduced by taking the expansion in a variable 
h = c X s for a suitable constant c, rather than in s = t^^'^'^. This trick was very useful 
for computing the larger solutions (with > 15 branches say). 
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4) To obtain the coefficients of the polynomial F from the Puiseux expansion we use 
the trick suggested in ^Oj: Write F in the form 

F = q{t)y'' + + ■ ■ ■ + pr{t)y + po{t) 

for polynomials pj, g in t and define rational functions rj(t) := pi/ q for i = 0, . . . , — 1. If 
yi, . . . ,yN denote the (locally defined) solutions on the branches then for each t we have 
that yi(t), . . . , yNif) are the roots of F{t, y) = and it follows that 

y"" + rM-i{t)y''-^ + ■■■ + n{t)y + ro(t) = - yi{t)){y - y^it)) • • ■ - y^it)). 

Thus, expanding the product on the right, the rational functions rj are obtained as sym- 
metric polynomials in the yf. 

ro= i-l)^yi---yN, ■■■ ,rAr_i = -(yiH h2/iv). 

Since the are global rational functions, the Puiseux expansions of the yi give the Lau- 
rent expansions at of the rj. Clearly only a finite number of terms of each Laurent 
expansion are required to determine each r^, and indeed it is simple to convert these 
truncated Laurent expansions into global rational functions. (This is easily done by Pade 
approximation, e.g. as implemented in the Maple command "convert ( ,ratpoly)".) 

5) Much time may be saved by carefully choosing the representative for the solution 
in the first place (i.e. try to choose an equivalent solution for which the polynomial F 
is simpler). Heuristically this can be estimated by seeing how complicated the algebraic 
numbers are (or by seeing how complicated the coefficients of the polynomial q{t) 
are; this is usually easily obtained from {yi + ■ ■ ■ + j/at), i.e. before having to compute 
complicated symmetric functions). 

6) Use Okamoto symmetries wherever possible: e.g. if (we can arrange that) the solution 
has the symmetry {y,t) i— > {y/t,l/t), swapping 6*2 and 6*3 then the coefficients of each 
Pi, q should be symmetrical, thereby essentially halving the number of coefficients that 
need to be computed. (Also for the 24-branch icosahedral solution in [S] it was too 
cumbersome to compute the longest symmetric functions, corresponding to the 'middle' 
polynomials Pi, but, by using another Okamoto symmetry, the outstanding coefficients 
could be determined in terms of those we were able to compute.) 

7) Finally there are various optimisations that can be made (especially in computing the 
symmetric functions of the Puiseux expansions) if we expect F to have integer coefficients 
(which is the case for all examples so far). 

Acknowledgments. The author is very grateful to Mark van Hoeij for help computing the more 
difficult parameterisations of the curves F = 0, and to both C. Doran and A. Kitaev for explaining 
various aspects of their work. 

After this work was complete A. Kitaev informed the author that he had found an explicit family of 
covers corresponding to the genus one (2, 3, 7) solution of section Eland had obtained a similar solution. 
Happily the solution here and that of Kitaev are not related by Okamoto transformations, but arise 
by choosing different embeddings of A237 into PSL2(C). In fact there are three inequivalent choices, 
corresponding to the three conjugacy classes of order 7 elements in PSL2(C) = S03(C). (This is analogous 
to the sibling solutions of 5 which arose from the two classes of order 5 elements.) The third inequivalent 
Pvi solution is: 

_ 1 (s^o + 5 + 24 + 20 - 266 - 2874 - 14812 - 40316 - 85359 - 100067 s - 67396) u 
^ ~ 2 16 (s + 1) (s^ + s + 7) (5 + 63 s5 + 252 s* + 854 + 1449 + 1827 s + 2030) 

with t,u, s exactly as in lO and = (4/7,4/7,4/7, 1/3). 
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